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ABSTRACT 

This paper analyzes astrophysical scenarios that may be detected at the 
upper end of the energy range of the Gamma Ray Large Area Space Tele- 
scope (GLAST), as a result of cosmic-ray (CR) diffusion in the interstellar 
medium (ISM). Hadronic processes are considered as the source of 7-ray pho- 
tons from localized molecular enhancements nearby accelerators. Two particular 
cases are presented: a) the possibility of detecting spectral energy distributions 
(SEDs) with maxima above 1 GeV, which may be constrained by detection or 
non-detection at very-high energies (VHE) with observations by ground-based 
Cerenkov telescopes, and b) the possibility of detecting V-shaped, inverted spec- 
tra, due to confusion of a nearby (to the line of sight) arrangement of acceler- 
ator/target scenarios with different characteristic properties. We show that the 
finding of these signatures (in particular, a peak at the 1-100 GeV energy region) 
is indicative for an identification of the underlying mechanism producing the 7- 
rays that is realized by nature: which accelerator (age and relative position to 
the target cloud) and under which diffusion properties CR propagate. 

Subject headings: 7-rays: theory, 7-rays: observations 
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1. Introduction 



1506. This is the number of cosmic photons with energies above 10 GeV that were 
detected during the 9 years lifetime of the Energetic Gamma Ray Experiment Telescope 
(EGRET, Thompson et al. 2005). Out of this number of photons, 187 photons were 
found within 1 degree of sources that are listed in the Third EGRET Catalog (Hartman 
et al. 1999) and can be plausible related to the more energetic extent of the cataloged 
EGRET sources. The majority of the remaining photons correspond to diffuse Galactic and 
extragalactic radiation, albeit this conclusion is based on the similarity of their spatial and 
energy distributions with the diffuse model. No significant time clustering nor source shining 
only at such high energies was detected. The scarcity of these detections represents the last 
electromagnetic window that remains to be opened between already explored energy ranges. 
EGRET was hampered in performing detailed studies of the 7-ray sky above 10 GeV, due 
to back-splash of secondary particles produced by high-energy 7-rays causing a self-veto in 
the monolithic anti-coincidence detector used to reject charged particles, and due to a non- 
calibrated detector response. GLAST, soon to be launched, will not be strongly affected by 
these effects since the anti-coincidence shield was designed in a segmented fashion (Moiseev 
et al. 2007). The effective area of the GLAST will be roughly an order of magnitude larger 
than that of EGRET leading to an increased sensitivity for detecting celestial 7-ray photons 
(see Fig. 1 of Funk et al. 2008). GLAST's default observing plan is a survey mode where 
the sky is uniformly observed. The increase of sensitivity and the survey mode open the 
gate to the possible discovery of new phenomenology. This paper analyzes astrophysical 
scenarios that can be detected at the upper end of the energy range of GLAST, as a result 
of cosmic-ray (CR) diffusion in the ISM. 



2. The role of diffusion 

The 7r°-decay 7-ray flux from a source of proton-density n p is given by 

POO 

F{E,) = 2 / (F^Ej/y/E* - ml) dE w , (1) 
where E™ m is the minimum pion energy given by E^ m (E 1 ) = E^ + rn^/AE^, and 

.glial 

F w (E n ) = 4vm p / " J p {E){da v {E v , E p )/dE w ) dE p , (2) 
Je™™ 

with da 7T (E 7T , E p )/dE n being the differential cross-section for the production of ir° of energy 
E n by a proton of energy E p in a pp collision. For an study of different parameterizations 
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of this cross section see, e.g., (e.g., Domingo-Santamaria & Torres 2005, Kelner et al. 2006). 
The limits of integration in the last expression are obtained by kinematic considerations 
(e.g. Torres 2004). We have implicitly neglected any possible gradient of CR or gas number 
density in the target. The CR spectrum is given by J P (E, r, t) = [c/3/47r]/, where f(E, r, t) 
is the differential number density of protons at an instant t and distance r from the source. 
Just for comparison, the spectrum of cosmic-rays in the Earth vicinity is (e.g., Dermer 1986) 
J P (E) = 2.2(£/GeV)- 2 - 75 cm-V 1 GeV- 1 sr- 1 . 

The particle's differential number density / satisfies the radial-temporal-energy depen- 
dent diffusion equation 

(df/dt) = (D(E)/r 2 )(d/dr)r 2 (df/dr) + (d/dE) (Pf) + Q, (3) 

where P = —dE/dt is the energy loss rate of the particles, Q = Q(E, r, t) is the source 
function, and D(E) is the diffusion coefficient, for which we assume a dependence on the 
particle's energy. The energy loss rate are due to ionization and nuclear interactions, which 
timescale is t pp , with the latter dominating over the former for energies larger than 1 GeV. 
The nuclear loss rate is P nuc = E/t pp , with r pp = (n p c no-pp)^ 1 being the timescale for the 
corresponding nuclear loss, k ~ 0.45 being the inelasticity of the interaction, and a pp being 
the cross section. 

Aharonian & Atoyan (1996) presented a solution for the diffusion equation with an ar- 
bitrary diffusion coefficient, and impulsive injection spectrum f in ^(E), such that Q(E,r,t) = 
N f in j(E)SfS(t). For the particular case in which D(E) oc E 5 and / in j oc E~ a , it reads 

f(E, r, t) ~ (N E- a /7cV 2 Rl i{ ) exp [-(a - l)t/r pp - (R/R di{ ) 2 } , (4) 

where 

R Aii = 2(D(E)t[exp(t5/r pp ) - l]/[M/r PP ]) 1/2 (5) 

stands for the radius of the sphere up to which the particles of energy E have time to propa- 
gate after their injection. In case of continuous injection of accelerated particles, Q(E, t) = 
Q E~ a T(t), the previous solution needs to be convolved with the function T(t—t') in the time 
interval < t' < t. If the source is described by a Heavside function, T(t) = Q(t), and for 
times t less than the energy loss time, f(E, r, t) = {Q$E~ a / AixD(E)r){2/ ^Tk) j^^^e'^dx, 
(Atoyan et al. 1995). We will assume typical values, a = 2.2 and 5 = 0.5. 

In the case of energy-dependent propagation of CRs, a large variety of 7-ray spectra is 
expected (e.g., Aharonian & Atoyan 1996, Gabici & Aharonian 2007, Torres et al. 2008). 
Diffusion of CRs have also been explored as an explanation for the high energy observations 
of the Galactic Center (e.g., Hinton & Aharonian 2007). We have systematically studied, 
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numerically producing over 2000 E 2 F distributions, the dependences with the involved pa- 
rameters. Table [1] summarizes the results both for an impulsive and a continuous accelerator. 
The influence of the total energy injected by the accelerator as CRs (W p ), the age of the 
accelerator (t), the medium density in which the CRs propagate (n), the diffusion coefficient 
of the medium (given at 10 GeV, -Dio), the distance from the accelerator to the molecular 
cloud where the 7-rays are produced (R), the density of the cloud (nci), its mass (Mci), and 
radius (r d ), and the distance of such system to the observer (d) are described. The dominant 
dependences are related with the age of the accelerator and the diffusion coefficient. These 
parameters both impact onto the CR distribution. The faster the diffusion, the farther the 
target can be from the accelerator, and still be subject to a significantly enhanced CR spec- 
trum, see Fig. [IJ We use a source (cloud) parameterized in units of M 5 = Mci/10 5 M Q 
and dkpc — dj 1 kpc. transition, defined in the case of an impulsive accelerator, is the age for 
which the timescale for the corresponding nuclear loss becomes comparable to the age of the 
accelerator itself. D trans i tion is the value of the diffusion coefficient for which the SEDs stop 
displacing in energy keeping approximately the same flux, as inferred from Fig. [IJ 

Setting, as an example, reasonable parameters for the energy injected by the accelerator 
into cosmic-rays (e.g., W p = 5 x 10 49 erg for an impulsive source and L p = 5 x 10 37 erg/s for 
a continuous one) and for the interstellar medium density (e.g., n — 1 cm -3 ), we have found 
several scenarios for the appearance of hadronic maxima produced by diffusion. Examples 
are shown in Fig. [2} for the two types of accelerators. We have found that two kinds of 
peaks at this energy regime are possible: those that are not to be detected by an instrument 
with the sensitivity of EGRET or MAGIC, and those that are not to be detected by an 
instrument like H.E.S.S. or VERITAS (the latter are not shown in the examples in Fig. [2]). 
We also note that the impulsive accelerator produces more steeper maxima. We find that 
maxima in the SED, hadronically produced as an effect of diffusion of CRs, are possible and 
not uncommon at the high-energy end, where they produce 7-radiation at a level of flux 
detectable by the LAT. 

Fig. [3] shows, as contour plots, the energy at which the maximum of the SED is found 
for the cases of impulsive acceleration of cosmic rays, at different distances, ages of the 
accelerator, and diffusion coefficient. The energy-dependent propagation effects underlying 
our expectations for the diffusion of CRs in the studied parameter space are clearly depicted 
there. The diffusion radius, for t <C t pp , is Rdu(E) = 2a/ D(E)t, so that at a fixed age and 
distance, only particles of higher energy will be able to compensate a smaller D w , producing 
SED maxima at higher D-values. The smaller values of Dio that we study are expected 
in dense regions of ISM (e.g., Ormes et al. 1988, Torres et al. 2008). It is interesting to 
note that for many, albeit not for all, of the SEDs studied, the maximum in E 2 F space is 
found at energies beyond the energetic range of GLAST. On the other hand, we also note 
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that the use of Fig. [3] for the the interpretation of a GLAST observational discovery of 
a 1 - 100 GeV maximum, provides interesting clues about the nature of the astrophysical 
system that generates the 7-rays. First, we find these SEDs in cases where the scenario does 
not predict detectable emission at the EGRET sensitivity, so that they will represent new 
phenomenology. Second, we see from Fig. [3] that the range of accelerator-target separations 
and ages of accelerator that would produce such a 1-100 GeV maximum is rather limited 
(see in Fig. [3] the narrow contours for maxima at such energies), which would lead to a 
direct identification of the source, in case such system is found in the vicinity of the GLAST 
detection. 

We now focus on the case of one unresolved but composite GLAST source. We thus 
consider two separate accelerator-cloud complexes that are close to the line of sight such that 
GLAST observe them as a single source, within its PSF. This kind of scenarios would yield 
the observational signature of an inverted spectrum. Fig. H] shows four possible inverted 
spectra. The two figures in the top (bottom) part are generated by an impulsive (continuous) 
accelerator. The SED characterizing the oldest (youngest) accelerator is shown by dashed 
(dot-dashed) lines in each of the scenarios. Even if EGRET could have been able to weakly 
detect some the inverted source models we simulated, it could not conclusively relate it to 
such phenomenon due to its large low-energy PSF. The counterpart at higher energies is a 
bright source potentially detectable by ground-based telescopes. Due to continuous energy 
coverage, GLAST is a prime instrument to track this phenomenology. The right panel cases 
show particular examples in which the detection of the source by an instrument with the 
sensitivity of EGRET is not possible at all. The inverted spectrum is less deep in these 
scenarios. Less pronounced V-shaped spectra can be obtained with concomitantly lower 
fluxes at TeV energies. 



3. Solution validity and timescales 

Essentially, if we use the equations given above to compute 7-ray fluxes, we are assuming 
that there is no significant cosmic-ray gradient in the target (i.e., the gamma-ray emissivity 
is constant throughout the cloud). This assumption may be valid when the size of the cloud 
is less than the distance to the accelerator, and the diffusion coefficient inside and outside 
the cloud are not significantly different (or even if they are, the proton-proton timescale is 
larger than the time it takes for cosmic-ray to overtake the whole cloud). 

Here we did not parameterize on the mass of the cloud, but rather on the value of 
M 5 /dl pc , i.e., the mass of the cloud in units of 10 5 M divided by the distance squared in 
kpc, which we now call A for brevity. In order for the solutions to apply, then, two conditions 
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need to be satisfied, r pp > t and the size of the cloud r c i (or simply r below) has to be such 
r c i < R (see Table 1 to refresh the meaning of variables, and do not confuse r, the size 
of the cloud, with R, the separation between the cloud and the accelerator). The density 
of the cloud, presented in terms of A, is n c i = 10 6 Ad\ / 'r 3 cm -3 . The nuclear loss rate is 
t pp ~ 100r 3 c / (Ad%. ) yr. The second condition (i.e, r < R) would be immediately satisfied 
if, say, r = R/x, with x > 1. Then, an astrophysical scenario (not unique!) which fulfills 
also the first condition (r pp > t) too is found when B = a/ (100/ A) (R 3 / (x 3 t)) > d kpc . To 
be viable, of course dk pc < 8 kpc, what can be used to define a minimum value of x, given 
an specific model with fixed A, R, and t. A final check can be done by comparing that 
the density obtained by replacing the previous inequalities, n c i < (10 8 /t yr ) cm -3 is in the 
range found in molecular clouds (which for the ages considered, is always fulfilled). Table 2 
shows the configurations used in the different figures of our work, together with the minimum 
value of x such that r < R, t pp > t, and dk pc < 8 are all maintained and the solutions used 
applicable. 

It is also interesting to discuss the different timescales within the cloud (see Gabici et 
al. 2007). In the former work, it was shown that dynamical and advection cloud timescales 
do not play a relevant role in this problem. To estimate the degree of cloud penetration by 
cosmic ray it suffices to compare the loss and diffusion timescales. The former was derived 
above. The latter can be written as 

— S(!)(^)li)" 

where the diffusion coefficient inside the cloud has been parameterized as 

In here, £ < 1 would account for a possible suppression of the diffusion coefficient inside 
the cloud as compared with that of the environment. For the proton-proton timescale to 
be larger than the time it takes for cosmic-ray to overtake the whole cloud, following the 
notation above, and taking an average cloud magnetic field of 10 /xG, the following condition 
applies 

R Dip 0.29 / g/10^G \ 05 

xA^lOVs"^ £ \E/GeV ) ' U 

With the parameters summarized in Table 2, it can be shown that plenty of clouds exist such 
that the diffusion timescale is shorter than the energy loss one at all energies, and especially 
at the ones we focus at the higher end of GLAST, and so we can neglect this effect. Our 
values of A and typical galactic distances make in general for not so massive clouds. Our 
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values of R make in general for not so large clouds. In this situation, only for low diffusion 
coefficient (e.g., 10 26 cm 2 s _1 ) timescales could become comparable if £ < 1; but then again, 
no significant suppression is expected when the environment has already such a low D w . 

4. Concluding remarks 

Compton peaks (of which the first example could have been found already, Aharonian 
et al. 2006) are not the only way to generate a maximum in a SED in the range of 1- 
100 GeV. A large variety of parameters representing physical conditions in the vicinity of 
a CR accelerator could produce a rather similar effect. Distinguishing between these cases 
would require multiwavelength information, search for counterparts, and modeling. If such 
a maximum is interpreted hadronically, as a result of diffusion of CR in the ISM and their 
subsequent interaction with a nearby target, the results herein presented constrain, given 
the energy at which the maximum of the SED is reached, the characteristics of the putative 
accelerator, helping to the identification process. Indeed, one of the most distinguishing 
aspects of this study is the realization that these signatures (in particular, a peak at the 
1-100 GeV energy region) is indicative for an identification of the underlying mechanism 
producing the 7-rays that is realized by nature: which accelerator (age and relative position 
to the target cloud) and under which diffusion properties CR propagate, as it is exemplified 
in Figure 3. In a survey mode such as the one GLAST will perform, it might also be 
possible to observe rather unexpected, tell-tailing SEDs, like those V-shaped presented here, 
if observed with instruments having a limited PSF, predictably leaving many Galactic sources 
unresolved. Indeed, we finally remark about the V-shaped spectra presented here that we 
have focused in the situation where from the same place in the sky we have a double peak 
structure. If instead two nearby sources at different energies are indeed resolvable by GLAST, 
i.e., displaced already in the GLAST range, and both above the threshold for detection, we 
would have a clear case of morphology change and spatially dislocation with energy what 
would certainly make the study also interesting. As noted recently 0, a double peak could 
also be seen -coming from the same place in the sky- if we consider two populations of CRs 
interacting with the same cloud, for instance, the usual E~ 2 ' 75 Galactic bath of cosmic rays 
(producing a photon-spectrum peaking in the GeV regime) and the escaped CRs from the 
nearby source (producing a photon spectrum peaking in the TeV range). This possibility 
would, however, be applicable only to the case of very young accelerators, e.g., 2000 years 
or so, with clouds significantly separated from it, and thus it is a less general scenario than 



1 In a poster by Gabici, Aharonian and Casanova in the Gamma-2008 meeting, July 6-11, 2008, Heidelberg, 
and to our knowledge yet unpublished. 
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what we have explored in this paper as the mechanism for V-spectra production. 
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Table 1: Dependence of the SED (E 2 F vs. E) on the different parameters. Imp. (cont.) 
stands for the impulsive (continuous) accelerator case. Dependences upon cloud parameters 
such as density (nci), mass (Mci), and radius (r c i) are obvious and related. 



parameter symbol and meaning 



effect on the E 2 F distributions versus E 



Accelerator 



W p : total energy injected for the accelerator as CRs 
L v : energy injected per unit time 
increasing t: age of the accelerator 



imp.: overall scaling, small effects in the range if in the typical range 10 50 -10 51 erg 
cont.: overall scaling, small effects in the range if in the typical range 

10 37_ 1Q 38 erg / s 

imp.: peak displaces to smaller energies for a fixed distance, 

until t > ttransition, and the peak displaces to smaller fluxes 

cont.: peak displaces to smaller energies and larger fluxes, for a fixed distance 



Interstellar medium 



n: density 
increasing Dio: 

the diffusion coefficient of the medium (at 10 GeV) 



negligible effects in the typical range 0.5-10 cm , since t pp ^ t. 

for a fixed age: curve displaces to smaller energies until Dio > ^transition 
where peaks generated by clouds at large separation, R, displace up 
and peaks generated by clouds at smaller R displace down in the SED 

for a fixed distance: curve displaces to smaller energies until Dio > ^transition 
where peaks generated by older accelerators (larger t) displaces down 
and peaks generated by younger accelerators (smaller t) displaces up 
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Table 2: Configurations used in this work and validity of the CR solutions used. 
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Fig. 1. — SEDs generated by CR propagation in ISM with different properties. Fluxes 
correspond to a cloud with M^/d 2 kpc = 0.5. Curve for D w = 10 26 , 10 27 , and 10 28 cm 2 /s 
are shown with solid, dotted, and dashed lines respectively. Sensitivities of EGRET (red) 
and GLAST (blue), H.E.S.S. (magenta) (survey mode and pointed observations with typical 
integrations), and MAGIC (yellow), are shown for comparison (see Fig. 1 of Funk et al. 
2008 for details on sensitivities). 
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Fig. 2. — Examples of the model predictions for a hadronic maxima in the 1-100 GeV 
regime. The left panel shows the predictions for a cloud scaled at M 5 /d\ pc = 0.025, located 
at 20 pc from an accelerator of 10 4 yr, diffusing with D w = 10 27 cm 2 /s. The right panel 
curve shows the predictions for a cloud scaled at M$/d\ pc = 0.08 located at 10 pc from an 
accelerator of 10 3 yr, diffusing with D w = 10 28 cm 2 /s. Increasing the ratio M 5 /c? 2 pc , the 
curves move up maintaining all other features. 
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Fig. 3. — For each combination of age and accelerator-target separation, for which more 
than two thousand spectra where numerically produced, the energy of the maximum of such 
spectra are shown in a contour plot. The color of the different contours corresponds to the 
range of energy where the maximum is found according to the color bar on the top of the 
figures. From left to right, plots are created for the case of an impulsive source injecting 
protons in a medium with _D 10 =10 26 cm 2 /s, 10 27 cm 2 /s and 10 28 cm 2 /s. 
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Fig. 4. — The parameters for the plots are as follows, (top left) the dashed curve on the 
left: t = 4xl0 5 yr, R = 5 pc, M b /d 2 kpc = 0.01; the dashed curve on the right: t = 10 4 yr, 
R = 20 pc, M 5 /d 2 kpc = 0.1; (top right) the dashed curve on the left: t = 2xl0 6 yr, #=100 
pc, M 5 /dl pc = 3; the dashed curve on the right: t = 4xl0 3 yr, R = 15 pc,M 5 /dl pc = 0.1; 
(bottom left) the dashed curve on the left: t = 2xl0 6 yr, R=15 pc, M 5 /d 2 kpc = 0.004; the 
dashed curve on the right: t = 10 3 yr, R = 5 pc, M^jd\ pc = 1; ( bottom right) the dashed 
curve on the left: t = 2xl0 6 yr, R = 40 pc, M 5 /d\ pc = 0.017; the dashed curve on the right: 
t = 6xl0 4 yr, R = 30 pc, M 5 /d\ pc = 2.5. D w is set to 10 26 cm 2 /s. R is the accelerator-cloud 
separation. 



